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Abstract 



A method for computed tomography, includes, in order to scan an object to be examined with a conical 
beam originating from a focus and with a matrix-like detector array for detecting the beam, the focus is 
moved, relative to the object to be examined, on a focal path around a system axis. A detector array 
supplies an output data corresponding to received radiation. The output data is filtered, and the filtered 
output data is backprojected three-dimensionally in order to produce at least one slice of a layer of the, 
which has a layer thickness. The slice represents absorption values obtained from the output data of a voxel 
belonging to the layer for the radiation from the beam. Filtering is carried out at least also in the direction of 
the tangent to the focal path belonging to a respective focal position, and normalization is carried out for 
each voxel considered 
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Beschreibung 

[0001] Die Erfindung betrifft em Verfahren zur Bildrekonstruktion fur die Computertomographie, aufweisend die Ver- 
fahrensschritte 

5 

a) zur Abtastung eines Untersuchungsobjekts mit einem von einem Fokus ausgehenden konusformigen Strahlen- 
biindel und mit einem matrixarligen Detektorarray zum Detektieren des Strahlenbiindels wird der Fokus relativ zu 
dem Untersuchungsobjekt auf einer Fokusbahn um eine Systemachse bewegt, wobei das Detektorarray der empfan- 
genen Strahlung entsprechende Ausgangsdaten in Form von Strahlen in Facherstrahlgeometrie liefert, 
10 b) die Ausgangsdaten werden gefiltert, und 

c) die Ausgangsdaten werden anschiieBend zur Erzeugung wenigstens eines Schnittbildes einer eine Schichtdicke 
aufweisenden Schicht des Untersuchungsobjekts riickprojiziert, wobei das Schnittbild aus den Ausgangsdaten ge- 
wonnene Absorptions werte der zu der Schicht des Untersuchungsobjekts gehorigen Voxel fiir die Strahlung des 
Strahlenbiindels reprasentiert. 

15 

[0002] Es sind Verfahren dieser Art bekannt unler dem Begriff "Filtered Backprojektion" bekannt geworden, die je- 
doch als 3D- Verfahren, d. h. im Zusammenhang mit einem matrixartigen Detektorarray, keine in der Praxis als ausrei- 
chend erachtete Bildqualitat erbringen, weil wegen des Kegeiformigen Rontgen strahlenbiindels sogenannte "cone beam 
artifacts" auftreten. 

20 [0003] AuBerdem ist an diesen Verfahren nachteilig, dass redundante Daten, wie sie bei Spiralabtastungen mit gerin- 
gem Tischvorschub durch mehrmaliges Bestrahlen ein und desselben Voxels entstehen, nicht genutzt werden, was zur 
Folge hat, dass die dem Untersuchungsobjekt verabreichte Strahlendosis nur unvollstandig zur Bildgebung genutzt wird. 
[0004] Weiter gibt es "Qberlegungen im Zusammenhang 2D- Verfahren zur Bildrekonstruktion so vorzugehen, dass vor- 
laufige Bilder in groBer Zahl durch "Filtered Backprojektion" aus Ausgangsdaten berechnet werden, die aus zur Bildre- 

25 konstruktion an sich nicht ausreichenden Abschnitten der Fokusbahn stammen, wobei die vorlaufigen Bilder erst in ei- 
nem zweiten Schritt zu einem endgiiltigen Schnittbild reformatiert werden. Diese 2D- Verfahren sind fur Detektorarrays 
mit groBer Breite, d. h. groBer Erstreckung in Richtung der Systemachse, weniger brauchbar, weil dann eine extrem 
groBe Anzahl von vorlaufigen Bildern zu verarbeiten ist, was selbst bei groBer zur \ferfugung stehender Rechenleistung 
problematisch ist. 

30 [0005] Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren der eingangs genannten Art anzugeben, das ein Stei- 
gerung der Bildqualitat ermoglicht. 

[0006] Nach der Erfindung wird diese Aufgabe geldst durch ein Verfahren mit den Merkmalen des Patentanspruchs 1. 
[0007] Demnach erfolgt die Filterung in Richtung der zu der jeweiligen Fokusposition gehorigen Tangente an die Fo- 
kusbahn. Es hat sich gezeigt, dass mit dieser Filterrichtung eine besonders hohe Bildqualitat erzielt werden kann, Der 

35 Wahi dieser Filterrichtung liegt die Erkenntnis zugrunde, dass das beschriebene, eine hohe Bildqualitat ermoglichende, 
auf vorlaufigen Bildern basierenden 2D- Verfahren sozusagen in ein 3D- Verfahren ubergehen wiirde, wenn die der Be- 
rechnung vorlaufiger Bilder zugrundeliegenden Abschnitte der Fokusbahn so stark verkurzt wiirden, dass sie nur noch 
eine einzige Projektion umfassen wiirden, deren Daten dann in Richtung der l&ngente an die Fokusbahn gefiltert wiirden, 
und dass erwartet werden konnte, dass ein solches 3D- Verfahren dann eine vergleichbar gute Bildqualitat wie das 2D- 

40 Verfahren ermogiichen wiirde. 

[0008] Das erfindungsgemaBe Verfahren lasst sich dann besonders leicht implementieren, wenn vor dem Filtern eine 
Umrechnung der in Facherstrahlgeometrie gewonnenen Ausgangsdaten in Form von Strahlen P(a, p, q) in in Parallel- 
strahlgeometrie vorliegende Paralleldaten in Form von Strahlen P(8, p, q) (azimutales "rebinning") bzw. P(8, p, q) (voll- 
standiges "rebinning", d. h. azimutales und radiales "rebinning") vorgenornmen wird. Dabei sind bezugnehmend auf Fig. 

45 3 

a der Fokuswinkel 
P der Facherwinkel 

q der der z-Koordinate entsprechende Zeilenindex des Detektorsystems, 
0 = a + p der Parallelfacherwinkel 
50 p = Rpsin(p) die dem Strahlabstand von der Drehachse (Systemachse) entsprechende Parallelkoordinate, und 
R F der Radius der Fokusbahn. 

[0009] Die Riickprojektion der Paralleldaten erfolgt nach einer bevorzugten Ausfuhrungsform der Erfindung in der 
Weise, dass im Zuge der Ruckprojektion fiir jedes Voxel (x, y, z) fiir jedes 0 e [0, k[ fur die Strahlen P(8 + krc, p, q) bzw. 
P(8 + kit , p, q), deren Projektion entlang der Systemachse durch (x, y) geht, die Summe 
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60 gebildet wird, wobei 

x, y, z die Koordinaten des jeweiligen Voxels sind, 

k eine ganze Zahl entsprechend der Anzahl der in die Rekonstruktion einbezogenen Halbumlaufe des Fokus, 
p die Parallelkoordinaten derjenigen Strahlen, deren Projektionen langs der Systemachse durch die Koordinaten (x, y) je- 
weiliges Voxels (x, y, z) verlaufen, 
65 p die Facherwinkel derjenigen Strahlen, deren Projektionen langs der Systemachse durch die Koordinaten (x, y) jewei- 
iiges Voxels (x, y, z) verlaufen, und 

h eine die Schichtdicke der in dem erzeugten Schnittbild dargesteilten Schicht des Untersuchungsobjekts bestimmende 
Gewichtungsfunktion ist, sowie 
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d eine Funktion ist, die gleich dem Abstand des jeweiligen Strahls von dem entsprechenden Voxel (x, y) ist bzw. von dera 
Abstand des jeweiligen Strahls von dem entsprechenden Voxel (x, y) abhangt. 
[0010] Die Schreibweise 



If 



bringt dabei zum Ausdruck, dass die Summenbildung wahlweise fur durch azimutales "rebinning" oder durch vollslan- 
diges "rebinning" gewonnene Strahlen erfolgen kann, wobei die FiLterung tangential zur Fokusbahn im Falle des azimu- 
talen "rebinnings" einer Filterung in p-Richtung und im Falle des vollstandigen "rebinnings" einer Filterung in p-Rich- 
tung ist. 

[0011] Infolge der Summierung sowohl iiber k als auch q ist dann sichergestellt, dass alle durch ein und dasseibe Voxel 
verlaufende Strahlen berucksichtigt werden und die dem Untersuchungsobjekt zugeruhrte Strahlendosis somit vollstan- 
dig genutzt wird. 

[0012] Eine besonders bevorzugte Ausfiihrungsform sieht vor, dass zur Ruckprojektion der Paralleldaten die auf die 
Sutnme H der Gewichte h 
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normierte Summe 
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gebildet wird. Diese Vorgehensweise ermoglicht eine nochmals verbesserte Bildqualitat, da eine mogliche t)berbetonung 
von Voxein, die von mehr Strahlen als andere Voxel getroffen werden, beseitigt wird und somit entsprechende Artefakte 
vermieden werden. Der CT-Wert des jeweiligen Voxels wird durch Summation uber 0 gewonnen. 
[0013] Das erfindungsgemaBe Verfahren kann beispielsweise Anwendung finden, wenn gemaB einer Variante der Er- 
findung die Fokusbahn eine Kreisbahn ist (Tomogrammabtastung). GemaB einer bevorzugten Variante der Erfindung ist 
die Fokusbahn jedoch eine Spiraibahn, die dadurch zustande kommt, dass der Fokus auf einer Kreisbahn tun die System- 
achse bewegt wird und gleichzeitig eine Relativbewegung zwischen Fokus und Untersuchungsobjekt in Richtung der 
Systemachse erfolgt. Auf Basis einer solchen Spiralabtastung konnen problemlos auch groBere \blumina des Untersu- 
chungsobjekts untersucht werden. 

[0014] Im Falle einer Tomogrammabtastung gilt fur k normalerweise k = 1 oder k = 2, im Falle einer Spiralabtastung 
wird k unter Berucksichtigung der pro Vollumlauf erfolgenden Relativverschiebung in Richtung der Systemachse so ge- 
wahlt, dass der abzubildende Bereich des Untersuchungsobjekts vollstandig erfasst ist. 

[0015] Die Erfindung wird nachfolgend anhand eines in den beigefugten schematischen Zeichnungen dargesteilten 
Ausfuhrungsbeispiels naher erlautert. Es zeigen: 

[0016] Fig. 1 in teils perspektivischer, teils blockschaltbildartiger Darstellung ein mehrere Zeilen von Detektorelemen- 
ten aufweisendes CT-Gerat, 

[0017] Fig. 2 einen Langsschnitt durch das Gerat gemaB Fig. 1, und 
[0018] Fig. 3 ein das "rebinning" veranschaulichendes Schaubild. 

[0019] In den Fig. 1 und 2 ist ein zur Durchfuhrung des erfindungsgemaBen Verfahrens geeignetes CT-Gerat der 3. Ge- 
neration dargestellt. Dessen insgesamt mit 1 bezeichnete Messanordnung weist eine insgesamt mit 2 bezeichnete Ront- 
genstrahlenquelle mit einer dieser vorgelagerten quellennahen Strahlenblende 3 (Fig. 2) und ein als flachenhaftes Array 
von mehreren Zeilen und Spalten von Detektorelementen - eines von diesen ist in Fig. 1 mit 4 bezeichnet - ausgebildetes 
Detektorsystem 5 mit einer diesem vorgelagerten detektornahen Strahlenblende 6 (Fig. 2) auf. In Fig. 1 sind der Uber- 
sichtlichkeit halber nur acht Zeilen von Detektorelementen 4 dargestellt, das Detektorsystem 5 weist jedoch, was in der 
Fig. 2 punktiert angedeutet ist, weitere Zeilen von Detektorelementen 4 auf. 

[0020] Die Rontgenstrahlenquelie 2 mit der Strahlenblende 3 einerseits und das Detektorsystem 5 mit der Strahlen- 
blende 6 andererseits sind in aus der Fig. 2 ersichtlicher Weise an einem Drehrahmen 7 einander derart gegenuberliegend 
angebracht, dass ein im Betrieb des CT-Gerats von der Rontgenstrahlenquelie 2 ausgehendes, durch die einstellbare 
Strahlenblende 3 eingeblendetes, pyramidenronniges Rontgenslrahlenbundei, dessen Randstrahlen mit 8 bezeichnet 
sind, auf das Detektorsystem 5 auftrifft. Dabei ist die Strahlenblende 6 dem mittels der Strahlenblende 3 eingestellten 
Querschnitt des Rdntgenstrahlenbundels entsprechend so eingestellt, dass nur derjenige Bereich des Detektorsys terns 5 
freigegeben ist, der von dem Rontgenstrahlenbiindel unmittelbar getroffen werden kann. Dies sind in dem in den Fig, 1 
und 2 veranschaulichten Betriebsmodus acht Zeilen von Detektorelementen 4, die im Folgenden als aktive Zeilen be- 
zeichnet werden. Die weiteren punktiert angedeuteten Zeilen sind von der Strahlenblende 6 abgedeckt und daher nicht 
aktiv. 

[0021] Jede Zeile von Detektorelementen 4 weist eine Anzahl K von Detektorelementen auf wobei = pi bis |5 K der 
Kanalindex ist und jedem Detektorelement ein Facherwinkel p k zugeordnete ist. Der Facherwinkel des mitderen Detek- 
torelementes ist gleich Null; die Facherwinkel der beiden au Bersten Detektorelemente sind pi = +p max und P& = -P ma x- 
[0022] Die aktiven Zeilen Lq von Detektorelementen 4 sind in Fig. 2 mit hi bis Lq bezeichnet, wobei q = 1 bis Q der 
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Zeilenindex ist der im Falle des beschriebenen Ausfuhrungsbeispiels sozusagen der z-Koordinate entspricht 

Kngswinkel des Rontgenstrahlenbiindels in einer die Systemachse Z und den Fokus F enthaltenden Ebene handelt 
SSSswinkeldesMntgenstrah^ 

enthaltenden Ebene (Facheroffnungswinkel) betragt2|W und ist in F.g. 1 eingetrage n . Svstemachse in R ot a- 

f 00241 Der Drehrahmen 7 kann mittels einer Antriebseinnchtung 22 urn eine mit Z bezeichnete Systemachse in Kota 
K^L^t werfen Die Systemachse Z verlauft parallel zu der z-Achse eines in Fig. 1 dargestellten raumhchen recht- 

iSS^ptrr D^ktorsystems 5 verlaufen ebenfalls in Richtung der z-Achse, wahrend^ ft , Z«ilen. deren 
Kb in Richtung der z-Achse gemessen wird und beispielsweise 1 mm betragt, quer zu der Systemachse Z bzw. der z- 

fo^l Ve Um2n Untersuchungsobjekt, z. B. einen Patienten, in den Strahlengang des RontgenstraUenbundelbringen izu 

LS^gt^^ 

KLS 'dSmtufintervalllon wenigstens » + 2 Pra „ (ein Halbumlauf plus F^ffnungswmkel) umfasst, das 
eine vollstandige Rekonstruktion eines CT-Bildes gestattet, wahrend em Vollundauf 2n umfasst 
•« ?00301 Die wahrend der Spiral- oder Tomogrammabtastung aus den Detektorelementen jeder aktiven Zei e des Detek 
SvstenSl pSefau geLenen, den einzelnen Projektionen P(a, |J, q) in FScherstrahlgeometne en^henden 
MesSwer^ 

do Serende Datenstrom zu einer Schnittbildrekonstruktionseinheit 13, die aus den Messdaten Schmttbilder von ge- 
^SSScSSS^ Untersuchungsobjekts nach einem noch im Einzelnen zu beschreibenden erfindungsgemaBen 

heSSen CverfUgung, die im Rahmen der Erfindung einer 3D-Ruckprojekdon "nterzogen werden. Ab End^b- 
50 nUstLnTD-BUddaten in Term einer dreidimensionalen Matrix beispielsweise mit den Achsen x, y, z zu Ve^gnng, 
wobet^^ 

wert enthalt Dieienigen Elemente der dreidimensionalen Matrix die den gleichen x-, y-, ^ 6 ""^ rt a™ 56 "', J"™ 
dann jeweUs panares Schnittbild der dem maBgeblichen x-, y-, oder z-Wert entsprechenden Schtcht des Untersu- 

^X^vo^ 

SS^S^SSSS^ ll£S^S^S^ wird von einer Generatoreinheit 17 mit den 
5£d£ SSSSsutai beispielsweise der Rohrenspannung U, v^o|t djj^d- nj 
wendigen Werte einstellen zu konnen, ist der Generatoreinheit 17 eine Steuereinheit 18 nut Tastatur 19 zugeordnet, die 

K was daS veranschaulicht ist, dass die Steuereinheit 18 mit dem Bildrechner 11 verbunden ist 
^ Hnter anderem kann die Anzahl Q der akdven Zeilen von Detektorelementen 4 und damit die Position der 
SnbVen^ ufd 6 Sgeste.lt werden? wozu die Steuereinheit 18 mit den Strahlenblenden 3 und 6 = ^ten 
V^™heUen 20 und 21 verbunden ist. Weiter kann die Rotationszeit x eingesteUt werden, die der Drehrahmen 7 fur 
rine SusSge Umdrehung benotigt, was dadurch veranschaulicht ist, dass die dem Drehrahmen 7 zugeordnete An- 

^"^^ Verfahxen auch in Facherstrahlgeometrie zu reali- 
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sieren, wird das beschriebene CT-Gerat vorzugsweise in einem Modus betrieben, in dem das erfindungsgemaBe Verfah- 
ren in Parallelstrahlgeometrie realisiert ist. 

[0038] Demnach werden die bei der Abtastung des fur die jeweilige Untersuchung relevanten Korperbereichs des Pa- 
tienten 8 durch Spiral- oder Tomogrammabtastung in Facherstrahlgeometrie gewonnenen Daten zunachst in an sich be- 
kannter Weise durch ein im Allgemeinen als "rebinning" bezeichnetes Verfahren in Daten in Parallelstrahlgeometrie urn- 5 
gewandelt. Diese Umwandlung beruht auf einer Umsortierung der in Facherstrahlgeometrie gewonnenen Daten derail, 
dass aus unterschiedlichen in Facherstrahlgeometrie aufgenoinmenen Projektionen Strahlen entnommen und zu einer 
Projektion in Parallelstrahlgeometrie zusammengefugt werden. In Parallelstrahlgeometrie reichen Daten aus einem In- 
tervall der der Lange TE aus, um ein vollstandiges Bild rekonstruieren zu konnen. Urn diese Daten gewinnen zu konnen, 
miissen nichtsdestotrotz zur Daten in Facherstrahlgeometrie aus einem Intervall der Lange n + 2P max zur Verfugung ste- to 
hen. 

[0039] Eine Projektion in Parallelstrahlgeometrie ist in Fig. 3 veranschaulicht. Demnach nehmen alle n ParaUelstrah- 
len RP t bis RPn dieser Projektion gegenuber der x-Achse des in Fig. 3 dargestellten und mit dem gemafi Fig. 1 iiberein- 
stimmenden Koordinatensystem den ParaUelfacherwinkel 0 ein. 

[0040] Anhand des in Fig. 3 in ausgezogener Linie dargestellten Parallelstrahls RPi soli im Folgenden der "Dbergang 15 
von Facherstrahl- auf Parallelstrahlgeometrie erlautert werden. 

[0041] Der Parallelstrahl RP X entstammt einer fiir die auf der Fokusbahn S liegende Fokusposition F t in Facherstrahi- 
geometrie gewonnenen Projektion. Der zu dieser Projektion in Facherstrahlgeometrie gehorige, durch die Drehachse 6 
und damit die z-Achse des Koordinatensystems verlaufende Zentralstrahi RFzi ist in Fig. 3 ebenfalls eingetragen. Der 
Fokusposition Ft entspricht der Fokuswinkel oti; dies ist der Winkel, den die x-Achse und der Zentralstrahi RFzi ein- 20 
schlieBen. Gegenuber dem Zentralstrahi "RFzi vveist der Strahl RP t den Facherwinkel P auf. Somit ist leicht erkennbar, 
dass fiir den ParaUelfacherwinkel 0 gilt: 

0 = <x + |J 



[0042] Der rechtwinklig zum jeweiligen Parallelstrahl gemessene Strahlabstand p von der Drehachse 6 bzw. der z- 
Achse ist durch 



p = RFsin(p) 



25 



30 



gegeben. 

[0043] Wie anhand des in Fig. 3 in verstarkter Linie dargestellten, durch die Drehachse 6 bzw. die x-Achse verlaufen- 
den Zentralstrahls RPZ deutlich wird, handelt es sich bei diesem Strahl um den Zentralstrahi einer in Fachergeometrie fiir 
die Fokusposition FZ unter dem Fokuswinkel otz aufgenommene Projektion in Facherstrahlgeometrie. Da fur den Zen- 
tralstrahi einer in Facherstrahlgeometrie aufgenominenen Projektion P = 0 gilt, wird deutlich, dass fur den Fall von Zen- 35 
tralstrahlen gilt: 

Je nachdem, ob ein azimutales oder vollstandiges "rebinning" durchgefuhrt wird, liegen die Parallelprojektionen in der 
Form 

P(<X,0,q) 40 
oder in der Form 
P(6, p, q) 

45 

vor, wobei 

a der Fokuswinkel 

(5 der Facherwinkel 

q der der z-Koordinate entsprechende Zeilenindex des Detektorsy stems, 

0 = a + 0 der ParaUelfacherwinkel 50 
p = Rpsin(p) die dem Strahlab stand von der Drehachse (Systemachse) entsprechende ParaUelkoordinate, und 
Rp der Radius der Fokusbahn sind. 

[0044] In einer einer ersten Ausfuhrungsform des erfindungsgemaBen Verfahrens entsprechenden, mitteis der Tastatur 
19 wahlbaren ersten Betriebsart arbeitet das beschriebene CT-Gerat auf Basis von durch azimutales "rebinning" gewon- 
nenen Projektionen. Die diescn Projektionen entsprechenden Daten werden im Falle des azimutalen "rebinnings" in p- 55 
Richtung gefiltert, d, h. jeweils in Richtung der zu der Fokusposition des Zentralstrahls der jeweiligen Parallelprojektion 
gehorigen Tangente T (siehe Fig. 3), und zwar unter Verwendung eines der in der Cornputertomographie iiblichen Filter- 
kerne, z. B. Shepp-Logan- oder Rainachandran-Lakshminarayanan-Kern. 

[0045] Die derart gefilterten Paralleldaten werden dann in der in der Weise ruckprojiziert, dass im Zuge der Ruckpro- 
jektion fur jedes Voxel (x, y, z) fiir jedes 0 e [0, n[ fiir die Strahlen P(0 + krc, p, q), deren Projektion endang der System- 60 
achse durch (x, y) geht, die Summe 

k it 

gebildet wird, wobei 

x, y, z die Koordinaten des jeweiligen Voxels sind, 

k eine ganze Zahl entsprechend der Anzahl der in die Rekonstruktion einbezogenen Halbumlaufe des Fokus, 
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P die Facherwinkel derjenigen Strahlen sind, deren Projektionen langs der Systemachse durch die Koordinaten (x, y) je- 
weiliges Voxels (x, y, z) verlaufen, und 

h eine die Schichtdicke der in dem erzeugten SchnittbUd dargestellten Schicht des Untersuchungsobjekts bestimmende 
Gewichtungsfunktion ist, sowie 
5 d eine Funktion ist, die gleich dem Abstand des jeweiligen Strahls von dem entsprechenden Voxel (x, y) ist bzw. von dem 
Abstand des jeweiligen Strahls von dem entsprechenden Voxel (x, y) abhangt. 

[0046] Infolge der gewahlten Filterrichtung und der infolge der Summation sowohl iiber k als auch uber q werden zum 
einen cone beam artifacts vermieden und zum anderen im Interesse einer hohen Dosisnutzung aller durch ein \bxel (x, y, 
z) verlaufenden Strahlen beriicksichtigt. 
to [0047] Der einem Voxel x, y, z zugeordnete Absorptionswert u x>y , z wird durch Summation uber 8 iiber mindestens ei- 
nen Halbumlauf gewonnen, d. h. durch Bildung von 

e 

[0048] Die dem jeweiligen Absorptionswert entsprechende CT-Zahl wird in herkommlicher Weise aus dem Absorpti- 
onswert ermittelt. 

[0049] Dabei konnen mittels der Tastatur 19 unterschiedliche Gewichtungsfunktionen h und unterschiedhche Funktio- 
nen d eingestellt werden . 

20 [0050] Als Gewichtungsfunktion h eignet sich beispielsweise eine Dreiecks- oder Trapezfunktion. 

[0051] Als Funktion d kann der Abstand des jeweiligen Parallelstrahls vom Voxel x, y, z oder stattdessen beispiels- 
weise die z(axiale)-Komponente dieses Abstandes eingestellt werden. 

[0052] In einer Abwandlung der vorstehend beschriebenen ersten Betriebsart wird zur Ruckprojektion der Parallelda- 
ten die auf die Summe H der Gewichte h 
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30 



normierte Summe 



H k q 

gebildet wird. Dies ermoglicht eine nochmals verbesserte Bildqualitat, da eine mogliche Uberbetonung von Voxeln, die 
35 in mehreren Halbumlaufen "beleuchtet", also von Strahlen getroffen werden, beseitigt wird und somit entsprechende Ar- 
tefakte vermieden werden. Diese Redundanz tritt bei Spiralabtastungen dann auf, wenn die pro Vollumlauf der Messan- 
ordnung stattfindende Relativverschiebung so gering ist (geringer Pitch), dass Voxel mehrfach bestrahlt werden. 
[0053] Eine einer weiteren Ausfuhrungsform des erflndungsgemaBen Verfahrens entsprechenden, mittels der Tastatur 
19 wahlbaren zweite Betriebsart unterscheidet sich von der ersten Betriebsart dadurch, dass das beschriebene CT-Gerat 
40 nicht auf Basis von durch azimutales, sondern von durch vollstandiges "rebinning" gewonnenen Projektionen arbeitet. 
Die diesen Projektionen entsprechenden Daten werden im Falle des vollstandigen "rebinnings" in p-Richtung gefiltert, 
also jeweils ebenfalls in Richtung der zu der Fokusposition des Zentralstrahls der jeweihgen ParaUelprojektion gehdri- 
gen Tangente T (siehe Fig. 3). 

[0054] Demnach wird fur die dcrart gefilterten Paralleldaten im Zuge der Ruckprojektion die Summe 



45 



gebildet, wobei 

50 p die Parallelkoordinaten derjenigen Strahlen sind, deren Projektionen langs der Systemachse durch die Koordinaten (x, 
y) jeweiliges Voxels (x, y, z) verlaufen. 

[0055] Auch im Falle der zweiten Betriebsart wird in einer Abwandlung zur Ruckprojektion der Paralleldaten eine auf 
die Summe H der Gewichte h 

normierte Summe, namlich die Summe 

tt k q 

gebildet. 

[0056] Im Falle der vorstehend beschriebenen ersten und zweiten Betriebsarten wird eine Funktionsweise des erfin- 
65 dungsgemaBen Verfahrens vorgesehen, bei dem im Zusammenhang mit einem Voxel (x, y, z) alle Strahlen betrachtet 
werden, deren Projektion entlang der Drehachse 6 bzw. der z-Achse durch x, y geht. Ob und in welchem MaBe diese 
Strahlen beriicksichtigt werden, wird durch die Gewichtungsfunktion h und die Funktion d bestimmt. 
[0057] Das CT-Gerat kann aber auch mittels der Tastatur 19 wahlbare weitere Betriebsarten aufweisen, die den zuvor 



6 



DE 101 59 927 A 1 

beschriebenen rait dem Unterschied entsprechen, dass fur eine gegebene Fokusposition der durch das jeweilige \bxel (x, 
y, z) verlaufende theoretische Strahl ermittelt wird und dann unter Beriicksichtigung der Gewichtungsfunktion h und der 
Funktion d in die Summenbildung im Zuge der Riickprojektion nur diejenigen Strahlen einbezogen werden, die tatsach- 
lich einen Beitrag zu der Summe liefem konnen, 

[0058] Im Falle der beschriebenen Ausfuhrungsbeispiele wird die Relativbewegung zwischen der Messeinheit 1 und 5 
Lagerungsvorrichtung 9 jeweils dadurch erzeugt, dass die Lagerungsvorrichtung 9 verschoben wird. Es besteht im Rah- 
men der Erfindung jedoch auch die Moglichkeit, die Lagerungsvorrichtung 9 ortsfest zu lassen und statt dessen die Mess- 
einheit 1 zu verschieben. AuBerdem besteht im Rahmen der Erfindung die Moglichkeit, die notwendige Relativbewe- 
gung durch Verschiebung sowohl der Messeinheit 1 als auch der Lagerungsvorrichtung 9 zu erzeugen. 

[0059] Das konusfbrmige Rontgenstrahlenbundel weist im Falle des beschriebenen Ausfuhrungsbeispiels einen recht- 10 
eckigen Querschnitt auf. Im Rahmen der Erfindung sind jedoch auch andere Querschnittsgeometrien moglich. 
[0060] Im Zusammenhang mit den vorstehend beschriebenen Ausfuhrungsbeispielen finden CT-Gerate der dritten Ge- 
neration Verwendung, d. h. die Rontgenstrahlenquelle und das Detektorsystem werden wahrend der Bilderzeugung ge- 
meinsam um die Systemachse verlagert. Die Erfindung kann aber auch im Zusammenhang mit CT-Geraten der vierten 
Generation, bei denen nur die Rontgenstrahlenquelle um die Systemachse verlagert wird und mit einem feststehenden 15 
Detektorring zusammenwirkt, Verwendung finden, sofern es sich bei dem Detektorsystem um ein mehrzeiliges Array 
von Detektorelementen handelt. 

[0061] Auch bei CT-Geraten der funften Generation, d. h. CT-Geraten, bei denen die Rontgenstrahlung nicht nur von 
einem Fokus, sondern von mehreren Foken einer oder mehrerer um die Systemachse verlagerter Rontgenstrahlenquellen 
ausgeht, kann das erfindungsgemaBe Verfahren Verwendung finden, sofern das Detektorsystem ein mehrzeiliges Array 20 
von Detektorelementen aufweist. 

[0062] Die im Zusammenhang mit den vorstehend beschriebenen Ausfuhrungsbeispielen verwendeten CT-Gerate wei- 
sen ein Detektorsystem mit nach Art einer orthogonalen Matrix angeordneten Detektorelementen auf. Die Erfindung 
kann aber auch im Zusammenhang mit CT-Geraten Verwendung finden, deren Detektorsystem in einer anderen Weise 
flachenhaftes Array angeordnete Detektorelemente aufweist. 25 
[0063] Die vorstehend beschriebenen Ausfuhrungsbeispiele betreffen die medizinische Anwendung des erfindungsge- 
maBen Verfahrens. Die Erfindung kann jedoch auch aufierhalb der Medizin, beispielsweise bei der Gepackpriifung oder 
bei der Materialuntersuchung, Anwendung finden. 

Patentanspriiche 30 

1. Verfahren zur Bildrekonstruktion fur die Computertomographie, aufweisend die \ferfahrensschritte: 

a) zur Abtastung eines Untersuchungsobjekts mit einem von einem Fokus ausgehenden konusformigen Strah- 
lenbiindel und mit einem matrixartigen Detektorarray zum Detektieren des Strahlenbiindels wird der Fokus re- 
lativ zu dem Untersuchungsobjekt auf einer Fokusbahn um eine Systemachse bewegt, wobei das Detektorarray 35 
der empfangenen Strahlung entsprechende Ausgangsdaten in Form von Strahlen in Facherstrahlgeometrie lie- 
fert, 

b) die Ausgangsdaten werden in Richtung der zu der jeweiligen Fokusposition gehbrigen Tangente an die Fo- 
kusbahn gefiltert, und 

c) die gefilterten Ausgangsdaten werden anschlieBend zur Erzeugung wenigstens eines Schnittbildes einer 40 
eine Schichtdicke aufweisenden Schicht des Untersuchungsobjekts dreidimensional ruckprojiziert, wobei das 
Schnittbild aus den Ausgangsdaten gewonnene Absorptionswerte der zu der Schicht des Untersuchungsob- 
jekts gehorigen Voxel ftir die Strahlung des Strahlenbiindels reprasentiert. 

2. Verfahren nach Anspruch 1, bei dem vor dem Filtern eine Umrechnung der in Facherstrahlgeometrie gewonne- 
nen Ausgangsdaten in Form von Strahlen P(a, P, q) in in Parallelstrahlgeometrie vorliegende Paralleldaten in Form 45 
von Strahlen P(8, P, q) bzw. P(0, p, q) vorgenommen wird, wobei 

a der Fokuswinkel 
P der Facherwinkel 

q der der z-Koordinate entsprechende Zeilenindex des Detektorsystems, 

8 = a + p der Parallelfacherwinkel 50 
p = Rpsin(p) die dem Strahlabstand von der Drehachse (Systemachse) entsprechende Parallelkoordinate, und 
R F der Radius der Fokusbahn sind. 

3. Verfahren nach Anspruch 2, bei dem die Riickprojektion der Paralleldaten in der Weise erfolgt, dass im Zuge der 
Riickprojektion fur jedes Voxel (x, y, z) fur jedes 0 e [0, n[ fur die Strahlen P(0 + krc, P, q) bzw. P(6 + kJT, p, q), de- 
ren Projektion entlang der Systemachse durch (x, y) geht, die Summe 55 



P 



q 



60 



gebildet wird, wobei 

x, y, z die Koordinaten des jeweiligen Voxels sind, 

k eine ganze Zahl entsprechend der Anzahl der in die Rekonstruktion einbezogenen Halbumlaufe des Fokus, 
p die Parallelkoordinaten derjenigen Strahlen sind, deren Projektionen langs der Systemachse durch die Koordina- 
ten (x, y) jeweiliges Voxels (x, y, z) verlaufen, 65 
P die Facherwinkel derjenigen Strahlen sind, deren Projektionen langs der Systemachse durch die Koordinaten (x, 
y) jeweiliges Voxels (x, y, z) verlaufen, und 

h eine die Schichtdicke der in dem erzeugten Schnittbild dargestellten Schicht des Untersuchungsobjekts bestim- 
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mende Gewichtungsfunktion ist, sowie 

deine Funktion ist, die gleich dem Abstanddes jeweiligen Strahls von dem entsprechenden \bxel (x, y) istbzw. von 
dem Abstand des jeweiligen Strahls von dem entsprechenden Voxel (x, y) abhangt. 

4. Verfahren nach Anspruch 3, bei welchem bci der Riickprojektion der Paralleldaten die auf die Summe H der Ge- 
wichte h 



f 



w 



normierte Summe 

H k q 



f 



9 + kn 



0 + kn, 



gebildet wird. 

5. Verfahren nach einem der Anspruche 1 bis 4, bei dem die Fokusbahn eine Kreisbahn ist. 

6. Verfahren nach einem der Anspruche 1 bis 4, bei dem die Fokusbahn eine Spiralbahn ist, die dadurch zustande 
kommt, dass der Fokus auf einer Kreisbahn um die Systemachse bewegt wird und gleichzeitig eine Relativbewe- 
gung zwischen Fokus und Untersuchungsobjekt in Richtung der Systemachse erfolgt. 
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